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Abstract. We show that the normahzation integral for the Schrodinger and Dirac scattering wave 
functions contains, besides the usual delta-function, a term proportional to the derivative of the 
phase shift. This term is of zero measure with respect to the integration over momentum variables 
and can be discarded in most cases. Yet it carries the full information on phase shifts and can 
be used for computation and manipulation of quantities which depend on phase shifts. In this 
paper we prove Levinson's theorem in a most general way which assumes only the completeness of 
states. In the case of a Dirac particle we obtain a new result valid for positive and negative energies 
separately. We also make a generalization of known results, for the phase shifts in the asymptotic 
limit of high energies, to the case of singular potentials. As an application we consider certain 
equations, which arise in a generalized interaction picture of quantum electrodynamics. Using 
the above mentioned results for the phase shifts we prove that any solution of these equations, 
which has a finite number of bound states, has a total charge zero. Furthermore, we show that in 
these equations the coupling constant is not a free parameter, but rather should be treated as an 
eigenvalue and hence must have a definite numerical value. 
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1 INTRODUCTION AND RESULTS 



It is a well known fact that the normalization integral of reduced radial Schrodinger wave 
functions contains a delta-function. A less well known fact is, that besides the delta-function, 
there are other terms. In this paper we shall show that the precise expression is 



dr Uki (r) Uk'i (r) = 2ti5 (k - k') + A {k, k') 

- (-1)' {2tx5 {k + k') cos [r]i (k) + r,i {k')] + A {k, -k')} , (1) 

where 

A(.,.') = {^';''w- , (2) 

rj'i (k) = dki]i (k) and rji {k) is a phase shift corresponding to the angular momentum / and 
linear momentum k. The corresponding result for the case of Dirac equation is 



dr [wieK (r) ui,'^ (r) + Ms^k (r) U2e'K {r)] = < 



27r6 {k -k') + A {k, k') 

- (-1)' {27r5 {k + k') cos [r/,« {k) + r/,« {k')] 
+A{k,-k')} , e = e' 
0, e^e' 

(3) 

where A{k,k') is the same function as in @, except that the Schrodinger phase shift rji 
is replaced by its relativistic counterpart r^^^- In most physical applications (for instance 
the elastic scattering cross section) the normalization integrals (0) or (j^) appear under the 
integral over k or /c', and hence the terms additional to the delta-function 6 {k — k') do not 
contribute. Yet these terms contain valuable information and can be used for the calculation 
of phase shifts or for the manipulation of quantities depending on them. In this paper we 
will use (P and (|^) to prove Levinson's theorem in a most direct and general way. 

This theorem is one of the most interesting nonperturbative results in quantum theory. 
It has many potential applications and has been applied recently in atomic physics 
in quantum field theories 0,0], and in solid state physics (where it is known in a modified 
form under the name Friedel's sum rule 0). In its original form the theorem says that 

r]i (0) = niTT , (4) 

which relates the scattering phase shift rji (0) at threshold (zero momentum) and for a given 
angular momentum / to the number of bound states ni of the Schrodinger equation with 
a spherically symmetric potential. If the Schrodinger equation has a zero-energy solution 
which vanishes at the origin and is finite at infinity and yet not normalizable (it is called 
a half-bound state or zero-energy resonance and is possible only if / = 0) then, as was first 
shown by R. Newton 0, Levinson's theorem is modified to read 

^o(0)= (no + ^)7r. (5) 

This result is subject to some restrictions on the potential. The most general proof of 
Levinson's theorem was carried out by Ni and Ma H] using the Sturm-Liouville theorem. 
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Their result for higher angular momenta / > 1 is 
resonance is possible, they obtained 



but for the S-state, where a zero-energy 



TT 



Vo (0) = noTT + ^ sin^ r]o (0) . 



(6) 



The authors did not try to solve this equation for rjo (0). However, we point out that this 
equation is easily solved and there are only three solutions: 



^0 (0) 




(7) 



The first solution is valid for the case without a zero-energy resonance, the other two if such 
a resonance exists. No examples are known for the second solution. In this paper we will give 
yet a more general derivation of the above results, which at the same time is considerably 
simpler. 

The first correct statement of Levinson's theorem for Dirac particles was given by Ma 
and Ni 

Vm. (0) + V-m. (0) = (iV+ + N-) 71 , (8) 

which is valid whenever there is no threshold resonance and 



(0) + v-n.. (0) = (iV+ + N-) n + i-iy ^ 



sin^ rimK. (0) + sin^ 7] 

—rriK (0) 



(9) 



which is valid for the case with a threshold resonance (which can appear only in the case 
K = ±1). Here ±m is the threshold energy of the Dirac particle, / = — 1 for k = — 1, —2, . . . 
and / = K for k = 1,2,... is the orbital angular momentum, is the number of positive 
and N~ the number of negative energy bound states of the Dirac equation and rj^rnn (0) are 
the phase shifts at threshold. Prior to the work of Ma and Ni claims were published stating 
that Levinson's theorem is valid for positive and negative energies separately and in the 
same sense as in the nonrelativistic case, i.e. ?7±mK (0) = N^tt, but later such claims were 
found incorrect |]10i- However, we shall prove in this paper that in a modified sense these 



claims are correct and that 



Vm,i (0) 



^0 + 4] ^' 

(< + 1) vr. 



V-mK (0) 




0,1, 





0,1, 







(10) 



'11^ 



where nf' and rij are the numbers of bound state solutions of certain Schrodinger equations 
which are given in the text and I = I — k/\k\. In (^) and (|TTp the first case refers to 
a situation without a threshold resonance and the other two cases to a situation with a 
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threshold resonance. Equations (|T^) and (|TT]) constitute the stronger statement of Levinson's 
theorem for Dirac particles. As a consequence of (||) and ([ToD , ([TlD it follows that 



n 



n-r 



(12) 



whereas in general A^^ 7^ ril and A*"" 7^ 11^ . 

Besides the phase shifts at threshold one can also obtain nonperturbatively the phase 
shifts in the asymptotic limit of high energies. A well known result is (see ref. [jll| p. 352) 



k—*oo 



— I dr V (r) 
k Jo 



where V (r) is a spherically symmetric interaction potential, and thus 



0, 



(13) 



(14) 



Obviously, equation (|TB|) is only valid if the integral on the right-hand side exists, which is 
not the case for potentials with a 1/r singularity or stronger. In this paper we shall extend 
this result to the case where the potential at the origin may be as singular as l/r^^° or less, 
including the obviously important case 1/r. At infinity the potential is assumed to vanish 
faster than 1/r. The result is 



dr 




(15) 



which still implies rji (00) = and which reduces to (|T3|) for potentials less singular than 1/r 
at the origin. In the case of a Dirac particle the Schrodinger result rji {00) = does not hold 
in general. Instead it was shown by Parzen that 



VeK {k) 



dr V (r) 



T-: I dr V (r) 



(16) 



where e = -^yW^+Tv? . In general, the right-hand side of ( |T6|) is a constant which yields 
?7±ooK (00) = only for a special type of potential. As was the case with (|13D the result ([T6|) 
only holds if the integral on the right-hand side exists. This is not the case for (r) ~ 1 /r at 
the origin. In this paper we shall derive (p!6| ) more rigorously than was done in |T^, and we 
shall generalize the result to the case of a potential behaving like 1/r at the origin. Unlike 
the Schrodinger case there is no reason to treat more singular potentials since the scattering 
wave functions do not behave well in this case. The result we shall prove is 



VeK {k) 



k 



dr 



'l-|T/(r)-l 



(17) 



which reduces to (|T^ if V is less singular than 1/r at the origin. 

As an application of the above results we shall investigate the following set of equations 



u 



K 



Uk + - Ulei^ - (mo + e + V?) U2e 



«2.K «2e« 

r 



- (mo 








(18) 
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+ = ^-,T.^ + . (19) 

0<eK<mo 0<eK<mo 

+ 7r~ (^) ^jEk (r') + Ui^-sK (r) Uj-sK (r')] , (20) 
where i = 1,2, j = 1,2, k = ±1, ±2, . . . , Ui^±s^n are bound state and Mj,±£K are scattering 



state solutions of (|IB|), e = ±e^, < < ""^O; are bound state and e = ±e, £ = yrriQ + k^, 
are scattering state energies and 

Qk= J2 + wL««) - E + "l-e.^) 

0<eK<"io 0<eK<mo 

f°° dk / 2 2 2 2 \ 

+ / + "2£K - - U2,-eJ ■ (21) 

For the sake of simphcity we assume that there are no threshold resonances and hence in (|20|) 
and ( pT] ) the region of integration excludes = 0. These equations are subject to certain 



boundary conditions which are explained in the text. Equation ([18|) is the radial Dirac equa- 
tion, (|TU|) is one of the Maxwell equations (Poisson equation) and ( PPD is the completeness 
relation for radial Dirac wave functions. Equations (p!8|-pO|) spherically symmetric 

special case of more general equations in a recently proposed generalized interaction pic- 
ture of quantum electrodynamics (QED). The derivation of these equations from QED goes 
beyond the framework of the present paper and we refer the reader to the forthcoming pa- 
per [0|. In the present paper, using the above results, we shall prove that for any solution 



of ( |T8[]20| ), which has a finite number of bound states, the total charge vanishes: 

Qo = 0, (22) 

where the charge density is defined by the right-hand side of (0) (in units of — Cq). Further- 
more, the coupling constant eQ/47r is not a free parameter but rather must have a numerical 
value for which 

dr (r) = . (23) 



2 NORMALIZATION OF SCATTERING STATES 
2.1 SCHRODINGER CASE 

Consider the radial Schrodinger equation 

^/(/ + 1) 



u 



kl 



+ 2mV - k'^ 



Mfci = , (24) 



for a scattering state characterized by the reduced radial wave function u^i (r) subject to the 
boundary conditions 

Mm(0) = 0, Mfc;(r) ^^2 sinUr-y + r/K^)) , (25) 
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where rji (k) is the phase shift. Here we assume that the potential V (r) is less singular at the 
origin than and that it vanishes at infinity faster than 1/r. The boundary conditions 
(l25|) determine the normalization of the wave functions: 



dr Uki (r) Uk'i (r) = 2'jt6 {k - k') + A {k, k') 

- (-1)' {2ti5 (k + k') cos [r,i (k) + r,i {k')] + A (fc, -k')} , (26) 

where 

A (fc, k') = I 2^'' ' 1 = I (27) 

and rj'i (k) = dkf]i (k). Notice that in the physical region where both k and k' are positive the 
last two terms in (^) vanish identically. Moreover, when integrated over positive values of k 
or k', only the first delta-function contributes. Thus in most cases one could drop all terms 
except the first. Then, however, one has lost the valuable information about the phase shifts 
Tji contained in the A {k, k') term. Therefore it is especially desirable to look into situations 
where only the diagonal terms in (p6|) are essential and the delta-functions do not contribute. 
For instance, subtracting from ( p6|) its noninteracting counterpart the delta-functions cancel 
and we obtain 



oo 

dr 



uli (r) - vl (r) = 2ri[ {k) + {-l)' 27r5 {k) sin^ r/, {k) , Vu (r) = 2krji {kr) , (28) 



where ji {kr) are the spherical Bessel functions. This equation turns out to be quite useful, 
as will be shown below. Also notice that the extra terms on the right-hand side of (|26|) 
cannot be cancelled by changing the normalization in (^3]). 

To derive equation (p6D we multiply (p^ ) by u^n and the corresponding equation for 
Uk'i by Uki- Subtracting the resulting equations from one another and integrating using 
Uki (0) = Uk'i (0) = 0, we obtain 

" ir u„ (r) u„ (r) = W W -"g W W , pg) 

where u'^^i {R) = dRUki {R). For large enough R one can evaluate the right-hand side of (|29|) 
using the asymptotic expression (^5|) and obtain 

1 

dr Uki (r) Uk'i {r) 

Using rji {—k) = —rji {k) (see ref. for a proof), equation (|26| ) follows as i? — cxo. 

It is instructive to consider a somewhat different derivation. Taking the derivative of (|2^) 
with respect to and multiplying the result by Uki, then multiplying (^) by dkUki, subtract- 
ing the resulting two equations from one another and integrating, we obtain 

J^"" dr uli (r) = ^ [u',i {R) dkUki {R) - um {R) Oku'^i {R)] . (31) 



2 Jo 
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For a large enough R we can use the asymptotics (^) to evaluate the right-hand side of (^Tj). 
The result ^ 

/ druli{r) = 2R + 2T^[{k)-{-l)^\sin[2kR + 2r]i{k)) (32) 
Jo k 

was first obtained by Liiders [|l^ . Equation (|32D is exact if the potential V (r) vanishes for 

r > R and is valid asymptotically otherwise. From (^) one easily recovers the diagonal 

terms of (|26|) . 



2.2 DIRAC CASE 

Consider now the reduced radial Dirac equations 



u'leK + - - {e + m-V) U2e 

r 



u. 



2eK 



r 



■U2eK + (e -m - V)ui^t, 





0, 



(33) 



where e is the energy, V (r) is the time component of a vector potential and k = ±1, ±2, . . . . 
The quantum number k is the standard parametrization of the total angular momentum 
j = |k| — 1/2 and of the relative orientation between the spin and the orbital angular 

ite boui 

UieK (r) 



momentum. The appropriate boundary conditions are 

Ml.. (0) = 

U2e. (0) = 



2 sin I /cr - y + r],^ (k) 



e + m Kok 2 Ur + rj,^ (k) 
2e e + m V 2 



(34) 
(35) 



where k = ±A/e^^-rr?, = I = \k\ — {1 — kq) /2, I = I — kq, and rj^^ {k) is the phase 

shift. To ensure the consistency of (0) and (|35| ) with (|33D we assume that V (r) behaves 
like or less singularly than 1/r at the origin and that it vanishes at infinity faster than 1/r. 
Thus the Yukawa potential e~^°^/r is allowed but the Coulomb potential is excluded. The 
normalization of the wave functions resulting from the boundary conditions ( ^4|) and ( ^5l) is 

'27T6{k-k') + A{k,k') 

r°°A I ^^ ^\ -{-ly {27r6 {k + k') COS [T],^{k) +T],^{k')] 

dr (r) Ml,,, (r) + Ms.k (r) Ms.'k (r)] = <^ (fc -k')} e = J 

[o, e^e'' 

(36) 

where A {k, k') is the same function as in (|27|), except that the phase shift rji is replaced by 
r]^f^. Notice the similarity of (^) to the corresponding Schrodinger case (^). As in that 
case we can extract the information on the phase shifts 77,, by subtracting from (|36D the 
corresponding equation for the noninteracting case, and we obtain 



dr 



u 



(r) + u 



2eK 



-'leK 



-'2tK 



2r^:jk) + {-iy2iTS{k)sm'7^,,{k) , (37) 



where 



le + m 



2e 



2krji {kr) , V2eK {r) 



It + m Kok 
2e e + m 



2kr jj (kr) 



(38) 
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is the regular free-particle solution of the Dirac equation (|33D , which obeys (|34D, (|35D and (|36|) 
with Tj^f^ ik) = 0. Contrary to the case of the Schrodinger equation the energy can now be 
either positive or negative and we can use instead of f to extract the diagonal terms 

in (EB). The result is 



(r) + mL (r) - ul_^^ (r) - ul _^^ (r) 

(k) - r/U. (k)] + i-iy 27r6 (k) sin [r/,. (k) + r/_,« (A;)] sin [r^,. (A;) - r/_,« (A;)] . (39) 



dr 



To derive ( pq ) we multiply the Dirac equation ( p3D by (— U2e'K, ^le's) and the Dirac equa- 
tion for e' by (— U2e/t, tiies) and subtract the results. Upon integration we obtain 

" dr (r) (r) + (r) (r)] = W W W W . (40) 

e — e 

For a large enough R the right-hand side of (^0]) can be evaluated using the asymptotic wave 
functions (EM) and (|35|). The result is 



R 



dr (r) (r) + MgeK (r) Mae'K (r)] 
sin [(A; - A;') R] 



(41) 



= A 

-(-l)'i? 
where 

A = 



k-k' 
r sin [(A; + k') R] 
1 VTk' 



r fi\ fi'\] , sin [r/e,. (A:) - ry,/^ jk')] 
cos [r7e«: (A;) - r],/^ (A; )J H ; cos [(A; - A; ) R\ 



cos [q^^ (k) + 1]^'^ {k')] + 



k-k' 
sin [r/eK (A:) + r],>^ {k')] 
k + k' 



cos [(A; + k') R] 



e' + m e 



tt' (e + m) 
B = 



e + m 



e + e' 



ee' (e' + m) / A; + A;' 



- — rki 



I e' + m 
ee' (e + m) 



e + m 



e + e' 



ee' (e' + m) Ik — k' 



(42) 



Equation ( ^61) is now easily recovered as i? ^ cxo. 

As in the Schrodinger case the diagonal terms of (pG]) can be obtained somewhat differ- 
ently. We derive the Dirac equation ( p3D with respect to e, multiply the resulting equation by 
(— M2eK, mick) and subtract from the result the Dirac equation multiplied by {—deU2eK, deUuK)- 
Upon integration we obtain 



dr 



ul^^ (r) + M2eK (r) = U2eK (R) C^eMles (R) " U^k (R) d^U2eK {R) ■ 



(43) 



Taking R large and inserting the asymptotic wave functions (^4|) and (^) on the right-hand 
side, we obtain 



R 



dr 



u 



leK 



(r) + u. 



2tK, 



I 771 

2R + 27]'^, {k) - (-1)' — sin [2kR + 2r/,« (A:)] 

6rL 



(44) 



From (44) one easily recovers the diagonal terms of (|36|) . Notice the close similarity of 
to the corresponding result for the Schrodinger case (P^. In particular for small energies 
(m/e 1) both results coincide. 
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3 HIGH ENERGY LIMIT 



3.1 SCHRODINGER CASE 

It is a well known result (see ref. [|ll[| p. 352) that 



m 



— I dr V (r) 



(45) 



and thus 

111 (oo) = . (46) 

Obviously, equation (^5]) is only valid if the integral on the right-hand side exists, which is 
not the case for potentials with a 1/r singularity or stronger. In this section we want to 
extend this result to the case where the potential at the origin may be as singular as 1/ r^"" 
or less, including the obviously important case 1/r. At infinity the potential is assumed to 
vanish faster than 1/r. The result we shall now prove is 



k—>oo 



k 



dr 



I 1m , , 



(47) 



which still implies (^Bf ) and which reduces to (^5)) for potentials less singular than 1/r at the 
origin. 

We introduce a new variable r and a new wave function x^/ (r) 



Uki (r) 



Xki \r) 



rr 

/ dr' F ir') 
Jo 



Fir) 



I 2m , , 



Putting this into the Schrodinger equation 

Xkl 



we get 



where 



A 



—2 9 

rv 9 



B 



V 



C 



V 



(4J 



(49) 



(50) 



and where the dots mean the derivatives with respect to the new variable r. Consider now 



a sphere of radius r = \Jm/k around the origin, k being large, and assume the case of the 
strongest allowed singularity V ~ l/r^"''. Then everywhere outside of the sphere the terms 
proportional to A, B and C in (^9]) are bounded by certain positive powers of m/k. Hence 
as — >■ cxD, we obtain 

+ 1) .2I 



Xkl 



(51) 



of which the physical solution is Xki {^) = 2krji (kr). Therefore as A; — 00, from ( ^ 
have 



we 



Uki [r 



17 ^^7 

2 sm \ kr \- k 



dr 



'1 



2m 



V{r)-1 



(52) 



and hence 
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3.2 DIRAC CASE 



In the case of a Dirac particle the Schrodinger result (^61) in general does not hold. Instead 
it was shown by Parzen [|12| that 



VeK {k) 



k 



dr V (r) 



— / dr V ir) 



(53) 



where e = itV/c^ + m?. In general, the right-hand side of ( p3D is a constant which yields 
^±ooK (oo) = only for a special type of potential. As was the case with (^), the result (|53|) 
only holds if the integral on the right-hand side exists. This is not the case for V {r) ~ 1/r 
at the origin. In this section we want to rederive (^) in a more rigorous way than was done 
in |T^ and to extend the result to the case of a potential behaving like 1/r at the origin. 



Unlike the Schrodinger case there is no reason to treat more singular potentials since the 
scattering wave functions do not behave well in this case. The result we shall prove is 



k^oo 



k dr 



(54) 



which reduces to ( ^3] ) if V is less singular than 1/r at the origin. 
Eliminating U2eK from (|33D, we obtain 



k(k + 1) K V 



Eliminating Mi^k from (^) leads to 



-V' + 2eV - k' 



UleK + 



v 



— Ml ,„ = , 

e + m-V 



(55) 



r e — m — V 



-V^ + 2eV - k^ 



U2eK + 



v 



e — m — V 



■ u 



2eK 



(56) 



Equations ( [55| ) and ( |5BD are equivalent to the Dirac equation (|33|) provided the boundary 
conditions (Q) and (|35| ) are imposed. Actually it is not necessary to solve (|55D and (|56| ) 
simultaneously. For instance, if (^) is solved for mi^^ one gets U2eK through the first of 
equations (|33|). In order to get rid of the last term in ([55| ) and ( p6D we introduce the new 
wave functions {wi^^, W2eK) defined through 



uien (r) 



'e + rri — y (r) 
2^ 



Wle^^ir) , U2eK{r) 



tkn e — m — V (r) 



\k\\K\ 



2e 



W2eK (r) . (57) 



Putting these in (|55D and (^), we obtain 



w 



la + i) K v 1 V" 

+ 



re + m — V 2e + m — V 



3 

+ - 



V 



4 \ e + m-V 



-V^ + 2eV - k^ 



wi,, = , (58) 
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W2.n = , (59) 



where we used k (k + 1) = / (/ + 1) and k (k — 1) = I (j[ + 1^ which are readily obtained 
from / = |k| — (1 — Kq) /2 and I = I — Kq, Kq = k,/\k,\. Equations (p^) and are of the 
Schrodinger-type except that the potential depends on the energy. Solving (^) and (|59|) is 
equivalent to solving the original Dirac equation (^) provided the boundary conditions 



Wu^ (0) = , WieK (r) 



2 sin (kr-^ + r]^^ {k) 



(60) 



W2en (0) = , W2eK (r) 



2 sin kr 



Til 



+ TjeK (k) 



(61) 



are met. We now turn to the limit k ^ oo. To perform this limit we can use either (|58|) 
or (pUD, the result will be the same. As in the Schrodinger case, considering a small sphere 
of radius (m/k)^^^, we observe that everywhere outside of the sphere the second, third and 
fourth terms in the parentheses of (|58D and (^) are bounded by positive powers of m/k as 
^ oo and hence can be neglected. The term V"^ behaves as 1/r^ at the origin and vanishes 
faster than 1/r^ at infinity. Hence it does alter the form of the wave function at the origin, 
but at large distances it contributes only to the phase of the wave function, which remains a 
Bessel function. However, V"^ is of order zero in k so that, according to the preceding section, 
it does not contribute to the asymptotic limit rj^^ (oo) and, in particular, cannot compete 
with the term 2eV which is of order one in k. Hence as k —>■ oo we neglect it and obtain 
from (HI) 

'/(/ + 1) 



+ 2eV - k' 



0, 



(62) 



which is subject to the boundary condition (^). Thus the high energy limit of the Dirac 
equation is dominated by the Schrodinger equation (24) in which the mass m is replaced 
by the energy e = ±VWTm?. It is easy to check that, despite this alteration of the 
Schrodinger equation, the arguments of the preceding section are still valid and we can 
^-Pply (0) I'splacing m by e. The result is given in (0). 



4 ZERO ENERGY LIMIT: LEVINSON'S THEOREM 



4.1 SCHRODINGER CASE 

The set of all physical solutions of the Schrodinger equation ( p4|) constitutes a complete set 
in the sense that for each I = 0,1, . . . 



POO (^J^ 

V Ueii (r) Ueii (r') + / —Uki (r) um (r') = 6 {r - r') 



(63) 
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where u^ii are bound state solutions of (0) normalized according to 



dr u 



(64) 



and Uki are scattering state solutions (|^). With respect to the completeness relation we 
have to distinguish two cases. In the first case, if there exists a zero-energy solution of the 
Schrodinger equation (^4]), which vanishes at the origin and is finite at large distances (we 
will call it finite), it is not normalizable in the sense of (|6^) . Such a solution, clearly, is not 
a bound state (it is called a half-bound state |Tl| or a zero-energy resonance since it 
becomes a bound state after an arbitrarily small increase in the strength of the potential). 
It is part of the continuum and is responsible for the delta-function in (pSj). Hence the sum 
over bound states on the left-hand side of ( |63D does not include the zero-energy solution and 
the integral over k is understood as including the value k = 0. This is the meaning of the 
completeness relation (^). In the second case, if a finite, zero-energy, solution exists, it is 
normalizable. Such a solution is a bound state and must be included in the sum over bound 
states on the left-hand side of (0) and, consequently, the integral over k is understood as 
excluding the value k = 0. The completeness relation for this case is 



E 



Ue,i (r) u,^i (r') + 



°° dk 
0+ 27r 



Uki (r) Uki (r') =5{r- r') 



(65) 



Notice that, if no finite zero-energy solution exists, then one can use either one of the 
completeness relations. In case a finite zero-energy solution exists, it is necessarily a zero- 
energy resonance if / = 0, and a bound state if / > 1. This is easily seen by examining the 
Schrodinger equation (p^) at large distances for potentials which vanish faster than 
and noticing that a finite zero-energy solution behaves as 1/r'. If the potential does not 
vanish faster than at infinity, then finite, zero-energy, solutions do not exist. 

Obviously the completeness relation also holds for a free particle, except that in this 
case there are no bound states and hence ( |63D and (|65D are identical. We now proceed with 
the second case above using (|65|) . Subtracting from (|^) the corresponding equation for the 
noninteracting case the delta-functions cancel and we obtain 



V Ueji (r) Ueji (r') + / — [uki (r) uh (r) - Vki (r) (r')] = , 
where Vki are the free particle solutions. Thus the diagonal part of (BH) reads 



(66) 



/•oo (j^j^ p 

E ^lii — uli (r) - vli (r) 



0. 



Integrating over r and substituting (|28|) , we obtain 

Vi (0) = Vi (oo) + niTT 



(67) 



(68) 



where ni is the number of bound states having angular momentum / (not counting the 2/ + 1 
degeneracy). Finally substituting (M), we get 



7]i (0) = niTT . 



(69) 
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Equation (pQ]) is known as Levinson's theorem |jl|]. It is valid whenever there is no zero- 
energy resonance and this, as explained above, is the case for all / > 1. This is, of course, a 
well known result. 

We now return to the first case where a finite zero-energy solution, if it exists, is not 
normalizable and where we have to use ( |63| ) instead of (]65|). Since a zero-energy resonance 
can occur only for / = we do not need to consider higher angular momentum states. Because 
the region of integration in (|63| ) contains the value k = 0, the steps preceding ( |69D remain 
the same except that now the delta-function on the right-hand side of (^) contributes, and 
we obtain 

Vo (0) = Vo (oo) + %7r + - sin^ r]o (0) . (70) 



Substituting (|46|), we get 



r]o{0) = no7T + -sm^r]o{0) . (71) 

This version of Levinson's theorem was first derived by Ni |0. Equation (^) can be solved 
for ?7o (0) and there are only three solutions: 



Vo{0) = { U + ijvr . (72) 




The second and third solutions correspond to the case with a zero-energy resonance, whereas 
the first solution corresponds to the case with no such state. The third solution is a result 
obtained by Newton 0. Equation (|72D is the most general statement of Levinson's theorem 
possible. The only input in our derivation is the completeness relation (^) and the boundary 
conditions (^). Although we do not have an example for the second solution in (|72|) there 
is little doubt that it can be realized. 

Levinson's theorem (|69[) can be proved in a number of different ways (see ref. 



and |TB[), for instance). The present proof has the advantage of being very general and 
extremely simple and also allows us some additional insights. 



4.2 DIRAC CASE 

A further advantage of the above proof is that it can be taken over to the case of the Dirac 
equation almost without change. In this section we will first rederive the result of Ma and 
Ni []TU|, which is the original correct statement of Levinson's theorem for Dirac particles. 



and which is valid for the sum of positive and negative energy phase shifts. Then we will 
prove a stronger statement of Levinson's theorem, valid for positive and negative energies 
separately. 

The zero-energy resonance of the preceding chapter is now called threshold resonance 
(its meaning is still the same: k = 0). As in the Schrodinger case we distinguish two cases: 
with and without threshold resonance. However, there is now a slight complication - the 
threshold resonance can have either a positive or negative energy, or there could exist two 
resonances simultaneously, one with positive and the other with negative energy. For the 
sake of simplicity we will consider only two cases; the gap will be filled at the end of this 
section where we consider the positive and negative energies separately. In the first case any 
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finite threshold solution (which does not necessarily exist) of the Dirac equation (|33D is a 
resonance and not a bound state. The completeness relation in this case reads 



— m<e„<m 
JO ZTT 



(73) 



where i = 1,2, j = 1,2, Uie^^ are bound state solutions of (^), £ = v^P~+"m2 and where 
the region of integration includes = 0, so that a term with a delta-function 6 (k) would 
contribute. The bound states are normalized according to 



dr 



«Lk (r) + (r) 



1. 



(74) 



In the second case any finite threshold solution (which does not necessarily exist) is a bound 
state. The corresponding completeness relation is 



6{r- r') 6ij = ^ Ui,^^ (r) uj,^^ (r') 

—m<e^<m, 

+ TT [uisK (r) UjeK {r') + Ui_e^ (r) Uj^^K (r')] , 



(75) 



where the region of integration does not include = 0, so that a term with a delta-function 
S {k) would not contribute. The rest of the derivation is identical to the Schrodinger case, 

). The result following from ( 



except that equation (|37D has to be used instead of 



nmn (0) + r^-™,. (0) = (iV+ + N-) vr , 



whereas the one following from ([73| ) is 

Vm. (0) + r/_„,. (0) = (iV+ + N-) 71 + (-1) 



I ^ 
2 



sin^ r]mK (0) + sin^ r] (0) 



IS 



(76) 



(77) 



where / = |k| — 1 for k = — 1, —2, . . . and / = k, for k = 1, 2, . . . is the orbital angular 
momentum, A*"^ is the number of positive and N~ is the number of negative energy bound 
states of (P^l). In deriving ( [77D we used ?7ook (oo) = — ?7_oo,k (oo) which follows from (Q). 
Equation ( [77D is Levinson's theorem for Dirac particles first obtained by Ma and Ni |[T0[| . 
Above we derived equations ( ^5] ) and ( pBD which, if subjected to the boundary conditions ( |5^ ) 
and (|35|), are equivalent to the Dirac equation (|33D. An inspection of these equations at 
threshold {k = 0) and for large r, shows that a threshold resonance can occur only if 
K = ±1. Therefore (|T6|) is valid for all k, = ±2, ±3, . . . and for k = ±1 if there is no threshold 
resonance, whereas ([77|) is valid for k = — 1 (/ = 0, I = 1) and k, = 1 (1 = 1,1 = 0). 

One may ask oneself whether (^6]) and ([77|) are perhaps valid for positive and negative 
energies separately. In fact, in some of the initial work (see refs. [|19| and 0) this was claimed 



to be true but later was found incorrect |]I0|. Yet, intuitively, by a number of reasons we 
expect that an extension of ( [76D and (ffTl), for positive and negative energies separately, should 
be possible. We now show that this is indeed the case. The basic observation which we need 
is the fact that the set of Schrodinger-like equations (|58| ) and (|59[) subject to the boundary 
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conditions (|60|) and ( |6TD is fully equivalent to the original Dirac equation (|33|) subject to 
the boundary conditions (|3^ ) and (^5|) . Equations (^) and (|59|) are useful because they 
are not coupled and hence the phase shift rj^^ (k) can be computed using any one of them, 
without reference to the other and for each of the positive and negative energies e separately. 
Certainly we cannot apply Levinson's theorem to (|58| ) or ( |59D directly since the potential in 
these equations depends on the energy and hence the completeness relation does not have 
the form (|63|). However, consider the following equations 



7/;+" - 
^kl 



1(1 + 1) K V 1 V" 3 

+ 7^7^ 77 + 



V 



r2m-V 21m -V 4V2m-F 



kl 



+ K V 1 V" ^3 



v 



r2m + V 22m + V 4:\2m + V^ 
which are subject to the boundary conditions 



+ 2mV - 



- 2mV - k^ 



w 



kl 



(0) = 



w 



kl 



2 sin kr 



(7J 



(79) 



(80) 



^:t(o) = o, 



w 



kl 

At threshold {k = 0) (^ and (|g 
and (ISTI) coincide with (Wi) and 



kl 



(r) 



2 sin (kr — ^ + (^) 



coincide with (|58D and ( |60D for e = m, and similarly ([79|) 
If) for e = — m. Moreover, both sets of equations and 



boundary conditions are analytical near the threshold. Therefore 

(0) = Vn.. (0) , v-7 (0) = r].^,. (0) . 



(82) 



Actually one would expect an ambiguity in both equations (|8^) , each in terms of an additive 
integer multiple of 27c. However, it is easy to see that both integers (say rii and 71,2) must be 
zero. This follows from the fact that the simultaneous change of m to —m and k to — k is 
a symmetry operation, which implies ni = n2. Equation ( |75D (or ([77[)), on the other hand, 
implies rii = —n2, and hence both integers are zero. Equations (|78|) and (ffOD are just usual 
Schrodinger equations, so that we can apply Levinson's theorem (|69|) and ( [72|) and obtain 



VmK (0) 




/ = 0,1, 

/ = 
/ = 



(83) 



(0) 



n-vr. 

Tin 



Tin 




0,1, 







(84) 



where n'^ is the number of bound state solutions (A;^ < 0) of (|78D and Uj is the number of 
bound state solutions of (|7Up. In ( P3| ) and ( |81D the first case refers to a situation without 
a threshold resonance and the other two cases to a situation with a threshold resonance. 
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Equations (|83D and (|8^ constitute Levinson's theorem for Dirac particles. As a consequence 
of and dH), (II we have 

Nt + N-=ni+r^, (85) 

whereas in general A^^ 7^ nf and N~ ^ ■ It is an interesting fact to notice that ([78| ) 
and (|79D do not correspond to the usual expansion based on the Foldy-Wouthuysen scheme 
(see [^). The above trick of "freezing" the energy of the second order Dirac equation was 



used for a different purpose in ref. 

5 APPLICATION TO A NONPERTURBATIVE QED 

Consider the following set of equations 

u'leK + - ^len - (mo + e + V?) U2eK = 

r 

«2e« ^2eK - (mQ - 6 - Ui^^ = , (86) 

r 

2 (=2 1 oo 

(5 (r - r') = ^ie.^ (r) Uje^K {r') + ^ -e.^ ('^) Uj-e^K {r') 

0<£K<mo 0<eK<mo 

+ TT (^) "^jEk {r') + Ui^^sK (r) Uj-en ijr')] , (88) 

J0+ 2,71 

where z = 1, 2, j = 1, 2, k = ±1, ±2, . . . , Ui^±e^n are bound state and Uj_±eK are scattering 



state solutions of (^61) , e = ie^, < £«; < ^o, are bound state and e = ±e, e = ymg + A;^, 
are scattering state energies and 



0<£re<»Tio 0<eK<mo 

POO (^f^ , 



o^ 



For the sake of simplicity we assume that there are no threshold resonances and hence in (|88|) 
and (P^D the region of integration excludes = 0. The boundary conditions are 



|¥?(0)|<oo, ¥?(oo)=0 (90) 
for the field if, and the usual boundary conditions 

Mi,±£«« (0) = , Ui^±e^^ (oo) — Ci^±e.K e"^"'' , (91) 



hold for the bound states, where Ci^±e^n are real constants and 7k = y — > g, and 



«i,K (0) = , Mi,K (r) W 2 sin (^r - ^ + r/,, (fc) j (92) 



16 



U2e. (0) = , W2e« (r) ^ 2 siii ^^^r - ^ + r/,. (fc)^ (93) 

for the scattering states, where / = \n\ — {1 — hq) /2,l = I — kq and kq = k/\k\. The bound 
states are normahzed according to 

«toW+4±..«W) = 1- (94) 



dr 

Notice that (|86| - p8|) is a closed system of equations. Also notice that any solution of this 



system constitutes a complete orthonormal set of functions and, therefore, Pauli's exclusion 
principle is built in. Equations (|86| - |88D arise as a spherically symmetric special case of 
more general equations in a recently proposed generalized interaction picture of quantum 
electrodynamics (QED). This generalized interaction picture is useful, since it allows for 
a new nonperturbative approach to QED. The derivation of these equations from QED 



and more details on their properties will be given in a forthcoming paper |13|. Here we 



want to apply the results of the preceding sections to the investigation of self-consistency of 



equations (86-p8D. We shall prove in this section that for any solution of (8^88]) with a finite 
number of bound states the total charge vanishes: Qo = 0, where the charge density is defined 
by the right-hand side of (0) (in units of —eg). Furthermore, the coupling constant el/iir 
is not a free parameter but rather must have a numerical value for which /q°° dr ip (r) = 0. 

Firstly we notice that there exists at least one solution. In fact, if we assume that there 
are no bound states then 

UUk = Vuk , U2eK = V2eK (95) 

is a solution of (BO-|SS|), where 



vie^ (r) = '^krji (kr) , v2eK (r) = \ j 2krjj{kr) (96) 

V 2e V 2e e + mo 



and ji (kr), jj{kr) are the spherical Bessel functions. In general, however, a solution of (^ 
will contain a number of positive and a number of negative bound states. While the 
possibility of an infinite number of bound states cannot be ruled out on simple grounds, we 
want to examine the consistency conditions for a finite number of bound states. 
If we multiply (^) by and integrate the resulting equation we obtain 

^'ir) = -f%-f-,f:^r dr' (rO + (r')] , (97) 

where 

Qo = -T.^ dr' [g^ (r') + g^, (r')] . (98) 

At large r we have 

0<£K<mo 0<eK<mo 
I <\k 1 

- (-1) 4m / cos [rje^ (k) - rj^sK ik)] cos [2A;r + rjet, [k) + rj^e^^ {k)] (99) 

J0+ ZTT e 
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and hence 



) ^2 I Jl ^2 I ^2 



+ (-1)' 2m / — — cos \r],^ (k) - r/_,, (A;)] sin [2kr + r/,« (A;) + r/_,, (A;)] . (100) 
J0+ Ztc ek 



T-^l^- + 0(e-) , a>0. (101) 



It now follows from (^) and (|100| ) that 

,2 

Now, if Qo > then ip acts as an attractive potential for the positive energy states and as 
a repulsive potential for the negative energy states. Similarly, if Qo < then tp acts as an 
attractive potential for the negative energy states and as a repulsive potential for the positive 
energy states. Consequently, in both cases ip supports an infinite number of bound states, 
like any attractive potential, which vanishes at infinity not faster than 1/r^. Thus Qo > 
and Qo < are inconsistent with a finite number of bound states. On the other hand, if 
Qo = then there will be at most a finite number of bound states, since ip acts as a short 
range potential with nonvanishing repulsive and attractive parts. The repulsive part is due 
to the electrostatic self-interaction of bound states, whilst the attractive part is provided 
by the scattering states, which fully compensate the repulsive part. Thus the consistency 
condition for the existence of a finite number of bound states is 

Qo = 0. (102) 

This consistency condition implies a restriction on ip. Combining (^Sj) and (P^D, we obtain 

°° { dk f°° r ^ 

(103) 

where N'^ and A^^ are the number of positive and negative energy bound states of ( |5B| ) 
respectively. Substituting (p9D, we get 

°° r 1 11 

Qo = - E l/^M^^-A^K -[^m«(0)-r/_^.(0)]- + [r/oo.(oo)-r7_oo.(oo)]- . (104) 

Since the boundary conditions for ip prohibit a singularity at the origin, we can use (^) for 
?7±oo (oo) and obtain 

Q^ = -J1 \^\[Nt-N--ni + T^ + -j^ dr^(r)|, (105) 

where we used Levinson's theorem (|8^), and where n/" and ii^ are the number of bound 
state solutions of ( [78D and (|79D with V = — ip and m = uiq respectively. Since the last term 
in (|105D is independent of k it must vanish, otherwise Qq would become infinite. Hence we 



conclude that 

/■oc 

dr (r) = (106) 
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must hold for any solution of (|8^-|88|) having a finite number of bound states. Notice that if all 



of the bound states have either positive or negative energies, then ( |106| ) implies (|102|) , which 
follows from (pS]). Also notice that /q°° dr ip (r) is a dimensionless quantity and therefore does 
not depend on the bare mass parameter mg and is a function of the coupling constant e^/Ai: 
only. Hence, due to ( |106| ), the coupling constant e^/An becomes an eigenvalue of ( |8^ - p8| ), 
i.e. we expect that ( |106|) holds only if e^/An acquires one or more specific numerical values. 
To illuminate this point we consider the following iteration procedure. As the initial step 
we make a reasonable guess of y?, which satisfies ( [L06| ) and which we denote by ipQ. Then we 
solve the Dirac equation ( PBD with (p = and obtain a complete set of bound and scattering 
states. At this stage c^/At^ is a free parameter. Now we fix it by requiring that with this 
complete set the solution of (0), which we denote by ipi, satisfies (|106|) . If there is no value 



of Cq/Att for which (|106|) is satisfied, then we start again trying to make a better choice of ipQ. 
If, however, there is such a value of Cq/Att, then we have completed the first iteration and 
can start the second one, now with pi instead of ipo. Assume that this iteration procedure 



is convergent, then it is clear that, as a result of (|106|) , the coupling constant will emerge, 
fixed at one or several numerical values. 
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